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Abstract 



We consider a burst of quadrupole gravitational radiation in the pres- 
ence of a large static mass M situated at its source. Some of the radiation 
is back-scattered off the static field of the large mass, forming a wave tail. 
After the burst, the tail is a pure incoming wave, carrying energy back 
] towards the source. We calculate this energy, and, in a numerical exam- 

ple, compare it with the outgoing wave energy. If M is sufficiently large 
the incoming energy can equal the outgoing energy, indicating that the 
primary outgoing wave is completely suppressed. 



1 Introduction 



It has been known for many years that gravitational waves have tails [1][2][3]. 
^ ■ (For a good literature survey see [4] ) . Tails are now recognised as factors in the 

planned observational detection of the waves [5] [14] [15]. The simplest sort of 
tail results from the back-scattering of an outgoing wave by the static field of 
the source, or of other bodies. Let us suppose for simplicity that there is an 
outgoing quadrupole wave that lasts for a finite time; when it ceases the tail 
persists, and as shown by Hunter and Rotenberg [3], it is then a pure incoming 
quadrupole wave. As such it carries energy, and this is what we calculate in 
this paper. The fate of this energy is not certain. It might be reabsorbed by 
the source, it might pass through the centre and become outgoing, or it might 
even form a singularity as suggested by the numerical work of Abraham and 
Evans [6] [7]. In this paper we do not consider these possibilities, but merely 
study the relative energies of the primary outgoing wave and the tail. 

The incoming energy carried by the tail must, of course, decrease the 
outgoing wave energy, as given, for example, by the quadrupole formula. The 
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question therefore arises whether the calculated incoming energy could be greater 
than the outgoing, which would presumably mean that the static fields of the 
source and its neighbours suppress the waves entirely From the theoretical point 
of view there seems no objection to this, since a similar suppression can occur 
of a light wave by a black hole. Indeed the suppression was investigated some 
years ago from a different point of view by Kundu [8] , and his work caused some 
debate [9] [10]. The conclusion seemed to be that a suppressive effect existed but 
would be counteracted by other phenomena, and would have no astrophysical 
importance [11]. 

The plan of the paper is as follows. We describe the metric and approx- 
imation method in Sections 2 and 3, and follow these with the linear approx- 
imation to the outgoing quadrupole wave in Section 4. We give the formulae 
for the wave tail in Section 5, and for the energy it carries after the vibration 
ceases in Section 6. Section 7 contains numerical calculations of the energies 
carried outward by the quadrupole wave and inward by its tail for a specified 
pulse. Conclusions are expressed in Section 8. 

2 Metric 

When dealing with outgoing waves we shall use the axially symmetric metric of 
Bondi [12]: 

ds 2 = -r 2 (Bd6 2 + B- 1 sin 2 9d(f> 2 ) + Ddu 2 + 2Fdrdu + 2rGd6du, (1) 

where B,D,F,G are functions of r,6 and u. u is a retarded time coordinate 
and may be thought of as ordinary time t at the source, which is represented 
by a singularity at the origin. The Schwarzschild solution in the coordinates of 

(1) is 

ds 2 = -r 2 (d8 2 + sin 2 9d<j) 2 ) + (1 - 2Mr- x )du 2 + 2drdu. (2) 

We shall also be studying incoming waves and for these we use the 
Bondi metric in the form 

ds 2 = -r 2 (Bd9 2 + B- 1 sin 2 8d(f> 2 ) + Ddv 2 - 2Fdrdv + 2rGd8dv, (3) 

with B,D,F,G now functions of r, 9 and v, and v an advanced time coordinate. 
In this metric the Schwarzschild solution is 

ds 2 = -r 2 (d6 2 + sin 2 9d<t> 2 ) + (1 - 2M r - 1 )dv 2 - 2drdv. (4) 

In these metrics the ranges are 

r > 0, < 9 < 7r, < 4> < 27T, -co < u, v < oo. 
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3 The approximation method 



We use an adaptation of the double series method. We suppose that the so- 
lution we are looking for can be expanded in a doubly infinite power series in 
M, a central spherical mass, and a parameter Q characterising the primary^] 
quadrupole moment of the vibrating source. The actual quadrupolc moment at 
time t is Qh(t), and h is zero except in t± < t < ti. Therefore we envisage an 
expansion of the gik of the form 

oo oo 

(P S ) 

the gik being independent of M and Q. 

If we substitute (5) into the vacuum field equations 

R ik = (6) 

and consider the coefficient of M P Q 3 we find seven equations of the form 

*im(K) = nth, (7) 

where the left-hand side is linear in the %k (and their derivatives), and the 

(qr) 

right-hand side consists of terms non-linear in the gik and their derivatives, 
with q < p,r < s and q + r < p + s. This is called the (ps) approximation. We 
can introduce this notation into the metric coefficients of (1), for example: 

-322 = r 2 B = r 2 [\ + Y.H M P Q S B ] 

(00) 

where B is zero. In this paper we shall go in detail only as far as the (11) 
approximation; this will express the interaction of the central mass and the 
quadrupole oscillations, which describes the simplest wave tail. To calculate 
the energies transmitted by the quadrupole wave and its tail we shall need 
formulae from the (02) and (22) approximations respectively which are known 
from previous work, namely eqns (12) and (29). 

x By primary quadrupole moment is meant the one inserted at the linear approximation, 
generating waves against a Minkowski background, and entering the usual quadrupole formula. 
Quadrupole waves appear in higher approximations and one might regard them as arising 
from higher quadrupole moments, though we do not do this here. Hereafter we drop the 
prefix 'primary', and refer simply to the quadrupole moment. 
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4 The linear approximation 



We write down the solution, in the coordinates of (1), of the linearised form of 
(6) (i.e. (7) with \1/ = 0) corresponding to a outgoing quadrupole wave [1]. 

°B = -sin 2 6(r- v h + r- 3 h), (8) 

D = -P 2 {2r- 1 h + 2r- 2 h + r- 3 h), (9) 

(01) 

F = 0, (10) 

(01) 1 

G = -sin6»cos6»(-2r- 1 /i + 4r- 2 /i + 3r~ 3 /i), (11) 

where P 2 (cos 9) is the Legendre polynomial and an overdot means differentiation 
with respect to u. h(u) is zero except for u\ < u < 112, where Ui = ti, and 
represents the quadrupole oscillation. We write Qh for the quadrupole moment 
of the source. We shall assume h is of class C 3 for u > u\ so that we can perform 
the necessary calculations on the tail in Section 7. 

The energy carried outward by this oscillation, calculated by the usual 
quadrupole formula or as a loss of mass as in [1], is 

7! 
30 



Eout = ^ / h du. (12) 



5 The wave tail 

The tail comes from the interaction of the quadrupole wave with the Schwarzschild 
metric. These are respectively the (01) and the (10) approximations, given by 
(8)-(ll) and (2). Inserting them into the right-hand side of (7) we have seven 

differential equations for the 'gik, of which the solution is [3] [13]. ^ 

C B = l s in 2 9(2r- 3 h + 10r- 4 h + r- 1 ^^Ar- 2 ^-+Ar- 3 H), (13) 
8 or 1 or 

(1 JD = -^P 2 (r- 3 h + 2r- 4 h + 2r- 3 H), (14) 

, ^ 

F = 0, (15) 
(11) 1 f)FT 

G = -sm28(3r- 3 h + 8,r- 4 h-2r- 2 — + 6r- 3 H), (16) 
8 or 

where H(r,u) is a tail function defined by 

H -~ 2 d-u]_ 00 (U + 2T-Q- (17) 

2 In the previous work cited the notation was different and the tail approximation was 
denoted (22) instead of (11). 
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From the definition of H we find 

dH dH n 

ir- 2 ir = 2r h ' 

or Ou 

now h = after the end of the quadrupole oscillation so 

H(r,u) = H(u + 2r), u > u 2 . (18) 

In the flat background spacetime u + 2r is the advanced time coordinate (here- 
after denoted by v) so for u > u 2 the tail junction H is an incoming wave. From 
(17) 

H _ h(u) g f u HQdt 

r J_ oo (u + 2r-0 2 

so as h as well as h is zero for ?t > « 2 , we have for the tail function after the 
vibration has finished 

* = 2 £<^f^r° > - <19 » 

Since Ho is a function of v = u + 2r we can write (13)-(16) for u > u 2 

as 

(1 B = isin 2 ^- 1 ^' -2r- 2 H Q + r- 3 H ), (20) 

D = -P 2 r- 3 H , (21) 
(ii) 

F = 0, (22) 

( G = ^sin29(-2r- 2 H + 3r- 3 H ), (23) 

where / means d/dv. This is the linear approximation for an incoming quadrupole 
wave in outgoing Bondi coordinates. It was given in [3], eqn (7.9) with s = 2 

2 

and /= —Ho- As v = u + 2r and u > u 2 , r > 0, (20)-(23) apply for 

v > u 2 , (24) 

u 2 being, of course, a constant. 

In order to calculate the energy carried by this incoming wave we write 
it in incoming Bondi coordinates, eqn (3). This too was done in [3] eqn (A25), 
and the result is 

°B = isin^r-^' + r- 3 ^), (25) 

( D = -P 2 (2r" 1 ^' -2r- 2 H[ ) +r- 3 H ), (26) 
(ii) 

F = 0, (27) 

( G = ^sm9cose(2r- 1 HH + Ar- 2 H -3r- 3 H Q ). (28) 
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6 Energy carried by the wave tail 



In the case of outgoing waves the energy carried is given by the loss of mass 
of the source. The simplest approximation in which it appears is that of the 
quadrupole-quadrupole interaction, so that in our notation of Section 3 it is 
of order Q 2 . This approximation was completely solved in [3], and the loss of 
mass, given by eqn (12), agreed with the ordinary quadrupole formula. 

The incoming waves described by (25)-(28) are of order MQ, but they 
have quadrupole form, and one can adapt the calculations of [3] to solve the 
M 2 Q 2 approximation in the incoming case. In fact the solution has exactly the 
same form as that of [3] 0, except for changes of sign. (See Appendix 1). The 
result is that the energy flowing inward after the vibration ceases is 

E in = ^f- J K 2 dv, (29) 

Hq being defined by (19). The lower limit of the integral arises from (24). 



7 An example 

To get a feeling for the magnitudes involved we construct an artificial example 
in which h(u) is a single pulse given by 

h(u) = (1 - u) 5 , -Ku<l; h(u) = 0, | u |> 1. (30) 

This has been chosen to make the tail function H reasonably simple. Although 
h is discontinuous at u = —1 it is very smooth at u = 1 which is sufficient to give 
an H which can be used in (29) if we put u% = 1. In this way we can calculate 
the energy carried inward by the tail after the end u = U2 of the vibration. 

First we calculate the energy E out carried outward by the pulse h ac- 
cording to (12): 

E out = % j h du = 768Q 2 . (31) 

where we have inserted the form (30) for h. 

From (19) we have for Ho(v), v > 1: 

H (v) = 2 / ) v > 1. (32) 

whence, after a fairly long but straightforward calculation, 

H'" = 16(l+i;)" 4 (64+30u-70i; 2 -90w 3 -30w 4 ) + 240(i;-l)[ln(w + l)-ln(u-l)]. 
(33) 

3 There described as the quadrupole-quadrupole approximation 
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To find Ei n we have to calculate the integral (29). This can be evaluated 
exactly using REDUCE. The result is 

E in = 28.2M 2 Q 2 (34) 

to three significant figures. 

These expressions for the inward and outward energies are in relativistic 
units. We need to find their ratio when they are expressed in units of customary 
dimensions. Denoting them in such units by Wi n and W ou t we show in Appendix 
2 that the ratio is 

W in _ G 2 m 2 H^dv 



Wout C 6 h w2 du ' 

Using the values of the integrals calculated above we find 

Wi 



(35) 



- 2 x lQ- 79 m 2 , (36) 



where m is in grams. To render Wi n comparable with W out the central mass 
would have to be that of a large black hole. 



8 Conclusion 

We have shown that in the presence of a sufficiently large mass the incoming 
energy of the tail can be comparable or even equal to that of the outgoing 
wave pulse. This raises the possibility that a pulse of radiation from, say a 
coalescing binary might be significantly reduced, or even suppressed, if the static 
gravitational fields in its neighbourhood are great enough. 

It must be emphasised that our calculations are not intended as a serious 
astrophysical model. We have treated the neighbours of the source as a single 
static spherical mass at the position of the source, and our numerical calculations 
are for one very special shape of pulse. Moreover, we have considered the tail 
only after the primary pulse has ceased. Nevertheless, we believe that the 
incoming tail energy should be taken carefully into account in gravitational 
wave experiments. 



9 Appendix 1 

In this Appendix we derive formulae relating to the incoming wave-tail energy 
referred to in Section 6. For this we make use of the outgoing quadrupole- 
quadrupole solution given by equations (3.7)-(3.11) in [3]. To obtain the com- 
plete M 2 Q 2 solution in the incoming case we simply make the coordinate trans- 
formation 

u = -v (37) 
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and the identification 

h = H (38) 

2 

with r, 9 and 4> remaining unchanged and where h is written as h in [3]. The 
only part of the transformation that is not entirely straightforward is the way 

in which Y of [3] transforms to its equivalent Y here. The former is defined 

by 

(24) 

where and overdot denotes differentiation with respect to u, and so Y must 
satisfy 

-§-ar <*>) 

where denotes differentiation with respect to v. Integrating with respect to v 
gives 

/>OC 

Y = - H'^dv. (41) 

J u 2 

The transformation (|37|), used on (3.7)-(3.11) of [3] then gives: 

2 35 A\ tjI ttIII \ ( o J f>„2 15 „4\ TTia-\-1 



r- 1 



+ [(-2 + 10s 2 - f s 4 )iT #o" + (-2 + 9s 2 - f s 4 )F ' 2 ]r" 

1/9 47 „2 , 331.4 m/ rr«„-3 , |Y 3 „2 27 4\tttt" 

+ + ^ 2 - |PK>- 4 + (|- 6s 2 + §* 4 )iWr- 5 

+(i - 3, 2 + f a *)flgr-« (42) 

( F = s *l-±H>*r-*-l 2 H HZr-*-l 2 H*r- 6 ] (43) 

(22) 

G = [(-^c S +§c S 3 )ilX V + (-!^ + Ic S 3 )^<-(^ CS +^ CS 3 )r]r- 1 
+ [(-2ca + ic8 a )H' Hg' + (-2cs + f cs 3 )iff ]r" 2 

,/27^ 185.J1TT/ pr«„-3 , rl.Jrr rjH , / 15 „„ , 105 „„3\ rr/2i„-4 
+ ( — cs - T6" cs ) H H r + l8 CS H H + (~T CS + "S -05 /""O J r 

+(f cs - ^cs^HoH^r- 5 + (§cs - f cs 3 )^" 6 (44) 

(22) 

5 = [(-^ a + ^VW + (-f^ + ^ 4 )W' + (-^ a -iio* 4 )5lr- 1 

+ | S 4 ^ 2 r- 2 + (-§ S 2 + fs^H^'r- 3 + [^H.H*' + - ^)H^ 

+ (-¥^ + fi^K^r" 5 + (-|s 2 + f S 4 )ff 2 r- 6 (45) 

(22) 

C = [(ij S 2 -^ 4 )!« + (is 2 -!T^ 

+i S 4 J? ' 2 r- 2 + (f S 2 - fs^H^r- 3 + [^H^H** + (-|§ s 2 + ^)H^ 



14" "U"u 1 V 16 1 32 



+ (f S 2 - f| S 4 )ff if r- 5 + (W ~ fs^Htr-z (46) 
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where s = sin6> and c = cos 6. By a coordinate transformation equivalent to 
(9.13) in [1], it can be shown that the above (22) solution describes, after the 
vibration, energy flowing inward equal to 

M 2 G 2 f°° 

E m = ^f- / H^dv. (47) 



10 Appendix 2 

We consider here the ratio of energies transmitted inward by the tail and out- 
ward by the primary quadrupole oscillation discussed in Section 7. 

We now write the energies in units of customary dimensions (ucd) and 
denote them by W out and Wi n . From (12) 

CO 2 f U2 2 

W out = ^ £ H du. (48) 

To express W in we first note that the dimensions of H are (time) _1 x (dimen- 
sions of h), so in ucd H becomes 

2 r mck 

We must also remember that M must be replaced by Gm/c 2 , where m is the 
central mass in ucd. Using this we have from (29) 

w - - <49 » 

Eqn (35) follows from (37) and (38) if we put u\ — — 1, Ui = 1. 
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